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Abstract 

We develop here a relatively simple description of dark energy based on the dynamics of non- 
minimally coupled to gravity phantom scalar field which, in limit, corresponds to cosmological 
constant. The dark energy equation of state, obtained directly from the dynamics of the model, 
turns out to be an oscillatory function of the scale factor. This parameterisation is compared to 
other possible dark energy parameterisations, among them, the most popular one, linear in the 
scale factor. We use the Bayesian framework for model selection and make a comparison in the light 
of SN la, CMB shift parameter, BAO A parameter, observational H(z) and growth rate function 
data. We find that there is evidence to favour a parameterisation with oscillations over a priori 
assumed linear one. 
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I. INTRODUCTION 



The recent discovery of the acceleration of the Universe is one of the most significant 
discoveries over last decade [l|, [2|. Observations of distant supernovae type la ^ as 
well as cosmic microwave background (CMB) fluctuations [l, Hj and large scale structure 
(LSS) jsl indicate that the Universe is undergoing an accelerating phase of expansion. These 
observations suggest that the Universe is filled by dark energy of unknown form, violating 
the strong energy conditions px + 3px > or a dynamical equation governing gravity 
should be modified. A simple cosmological constant model of dark energy can serve the 
purpose of explanation of dark energy and is in good agreement with the astronomical data 
(supernovae type la and other measurements). Although this model is favoured by the 
Bayesian framework of model selection j^-^l, it faces the serious problem of fine tuning 
[lO|. Therefore the other alternatives [uj have been proposed which includes an evolving 
scalar field. When one tries to accommodate a time-varying equation of state, the simplest 
parameterisation is the one which adds a linear dependence on the scale factor a. Other 
choices are motivated by a possibility of integration dark energy density in an exact form. In 
this context a class of simple oscillating dark energy equation of state coefficients appeared 
It is interesting that these models may provide a way to unify the early inflation 



and the late time acceleration. Moreover in these scenarios we obtain a possible way to solve 
the cosmic coincidence problem [l7- lol . 

If we allow that the dark energy density might vary in time (or redshift) then there appears 
a problem of choosing or finding an appropriate form of parameterisation of the equation 
of a state parameter wx{z). In the most popular approach wx{z) = px/px appears to be 
virtual and the dynamical dark energy parameterisation makes the model phenomenological, 
containing free parameters and functions. As a result we have a model difficult to constrain 
"201 . Another approach is to postulate a quintessence potential of the scalar field which 



has motivation from fundamental physics (particle physics) and then extracts it from the 
true dynamics directly [2l|-24|. In this approach we can expect that the parameterisation 
of dark energy equation of state refiects some realistic underlying of the physical model. 
The most popular dynamical form of dark energy offers the idea of quintessence. In this 
conception dark energy is described in terms of the minimally coupled to gravity scalar 
field with the potential V{(f)). The scalar field is rolling down its potential starts to 
dominate the energy density of the standard matter [25, 2^. The oscillating scalar field as a 
quintessence model for dark energy has been recently proposed [27-2^. The case of extended 



quintessence introduced by Amendola [30] was also considered in our previous papers [21 



231 ] where we assumed the non-zero coupling constant. The possibility of violating of the 
weak energy condition (phantom scalar field) was admitted. In this scenario, instead of 
standard minimally coupled scalar field, the phantom scalar field, non-minimally coupled to 



gravity, causes the accelerating phase of expansion of the Universe [21|, |22j . We found that 



in the generic case trajectories are approaching the de Sitter state after an infinite number 
of damping oscillations around the mysterious wx = — 1 value. Therefore the ACDM model 
appears as a global attractor in the phase space ip') (where ip is a phantom scalar field 
and ' = d/d In a). 

In this letter, we aim at testing and selecting the viability of different parameterisations 
of oscillating dark energy in the light of recent astronomical data. We focus our attention 
on the equation of state for a non-minimally coupled to gravity phantom scalar field with 
the potential in the simple quadratic form. 
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II. CLASS OF KINESSENCE MODELS 



This class of models is understood as a class of FRW models with standard dust matter 
and dark energy parameterised by redshift, i.e. wx = wx{z) 26]. For simplicity the 
flat FRW model {k = 0) is assumed. Then dynamics of the model is determined by the 
acceleration equation 



a 
a 



(p^ + (1 + 3Wx)px) = -\h^ [^m + (1 + 3Wx)nx 



(1) 



where a is the scale factor, dot is a differentiation with respect to the cosmological time, VL^ 
and VLx are density parameters for matter and dark energy X, respectively, H = (In a)' is 
the Hubble parameter. 

We assume that standard matter with energy density pm = Pm,o'^~'^ is a dust mat- 
ter and energy density of the dark energy is given, from conservation condition, by 

Px = px,oa'^ exp - 3 '^^^^^ da' . 

The acceleration equation ([T]) admits the first integral (which is called Friedman's first 
integral) in the form 

H " 



{i + zy + nx,of{z) 



(2) 



where Hq and Qifl are parameters referring to the present epoch, and z is the redshift related 
to the scale factor by the relation 1 + z = (the present value of the scale factor oq = 1). 
The phenomenological properties of dark energy are described in terms of the function f{z) 

such that f{z) = exp 3 ^'^I'+j''' ^ ^2;' . 

In the context of the accelerated expansion of the Universe most theoretical models of 
dark energy are based on scalar fields. It is a consequence of exploring an analogy to 



the inflationary theory of the primordial universe [3l|. However, a single canonical scalar 
field cannot explain the range of the coefficient of the equation of state w < —1 which is 
preferred by the astronomical data [8j]. One possibility, that has received much attention, is 
that we formally allow the scalar field to have a negative kinetic energy and switch its sign in 
comparison with canonical scalar field. There are some physical motivation for introducing 
such a phantom scalar field arising from string/M theory and in supergravity 32[. Another 
possibility lies in introduction of the coupling term ^Rip'^ between the scalar field and the 
gravity (for the review and references see [33]). Such a theory offers scalar-tensor models of 
a dark energy called extended quintessence. 

If we assume that a source of gravity is the phantom scalar field ip with an arbitrary 
coupling constant ^ then the dynamics is governed by the action 



S = \j d'xV^(mlR+{g 



(3) 



where 



(SttG) ^ and V{ip) is a scalar field potential. 



The phantom cosmology with general potentials was studied by Faraoni [3J] using lan- 
guage of qualitative analysis of differential equation to obtain the late time attractors without 
the specific assumptions of the shape of potential functions. There are many reasons why 
we should consider a nonzero coupling constant ^. First, a nonzero ^ is generated by quan- 
tum corrections even if it is absent in the classical action (see 35[ and references therein). 
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FIG. 1: The phase portrait represents generic behaviour of the system @ around the critical point 
of a stable node type. 



Another reason is that the non-minimal couphng is motivated by the renormahzation of the 
Einstein-Klein-Gordon equation. Of course the value of the coupling constant should be 
fixed by the physics only, but in relativity any value of the parameter ^ different from 1/6 
(conformal coupling) gives rise to the violation of the equivalence principle j36|. 

In our paper [22| we studied generic features of the evolutional paths of the flat FRW 
model with the phantom scalar field non-minimally coupled to gravity by using dynamical 
systems methods. We reduced dynamics of the model to the simple case of autonomous 
dynamical system on the invariant submanifold [ip^ip') (a prime denotes differentiation with 
respect to the natural logarithm of the scale factor) 



ip' = y 



-y - y^{y + 6^^) 
1 + (1 - 60y' 



1 + - 60 



(4) 



l + 6^?/'2(l-60 



where \^(^/') cx tl?' is assumed, and we found that principally there is one asymptotic state, 
which corresponds to the critical point in the phase space t/'q = il/v^W and il)'^ = 0. This 
critical point is also the de Sitter state [wx = ~1)- Note that in this model a problem of the 
big rip singularity does not appear in contrast to the standard phantom cosmology where it 
is present because the late time attractors in the phase space represent the de Sitter stage. 
There are two types of evolutional scenarios leading to this Lambda state (depending on the 
value of 0; through 

1. the monotonic evolution toward the critical point of a node type for < ^ < 3/25, 
(Fig. d]), in the special case ^ = 3/25 we obtain a degenerate node; 

2. the damping oscillations around the critical point of a focus type for 3/25 < ^ < 1/3, 

(Fig. ED. 

The effect of a non-minimal coupling can be treated as an effect of fictitious fluid with 
some effective coefficient of the equation of state given by 

(5) 
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FIG. 2: The phase portrait represents the generic behaviour of the system ^ around a focus type 
critical point. 

In both evolutional scenarios we can find linearised solutions of the dynamical system in 
the vicinity of the critical point, which the following Hartman-Grobman theorem are good 
approximations of the system. 

Finally one can compute linearised formulas for w{z) around the corresponding critical 
point for both cases (see Appendix 

In what follows we concentrate on the special case of parameterisation flA2|) with ^ = 1/6 



which corresponds to the conformally coupled phantom scalar field [36[. We will confront it 
(using the Bayesian model selection method) with most popular dark energy parameterisa- 
tions of Wx{z), which are presented in Table [H 

Recently, it has been argued that models with oscillating dark energy are favoured over 



a model with a linear parameterisation of EoS [17|, |23 



III. RESULTS 

A. Bayesian method of model comparison 

To find the best parameterisation of wx{z) we use the Bayesian method of model com- 



parison [39|. Here the best model (M) from the set of models under consideration is the 
one which has the greatest value of the probability in the light of the data {D) (posterior 
probability) 

, , , P(D\M)P(M) 
P{M\D) = ^ ^^^^^ ^ (6) 

P{M) is the prior probability for model M, P{D) is the normalisation constant and 
P{D\M) is the model likelihood (also called evidence) and is given by P{D\M) = 
f P{D\9^M)P{e\M)d9, where P{D\9,M) = L{9) is the likelihood function for model M 
and P{9\M) is the prior probability for the model parameters 9. It is convenient to 
consider the ratio of the posterior probabilities for models which we want to compare 
p|*j^j^| = p(j)[a^^) p(M2) " have no prior information to favour one model over an- 

other one (P(Mi) = P(M2)), posterior ratio is reduced to the ratio of the model likelihoods, 
so called Bayes factor (-B12), which values can be interpreted as the strength of evidence to 



favour model Mi over model M2 40[: < lni?i2 < 1 - 'inconclusive'; 1 < \nB12 < 2.5 
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TABLE I: Different dark energy parameterisations in terms of wx{z) = px/px ^ the coefficient of 
EoS. 



case 



parameterisation 



Chevallier-Polarski-Linder [37, 381 

Wxjz) =Wo + Wij^ 



(2) 



purely oscillating dark energy 

a) wx{z) = wq cos (wc In (1 + z)) 

b) wx{z) = -1 + Wo sin (ug In (1 + z)) 



(3) 



damping osc. DE 

a) wx{z) = wo{l + z)'^ cos (wc ln(l + z)) 

b) wx{z) = -1 + wo{l + zf sm{u}sln{l + z)) 



(4) 



damping osc. DE parameterisation determined directly from 
the dynamics of phantom scalar field model [22j 

wx{z) = -1 - 1(1 + (^(cos (^ In (1 + z)) + ^ sin (^ In (1 + z)))xo+ 
(cos (^ In (1 + z)) + ^ sin (^ In (1 + z)))yo) 
-|(1 + z)3 (^(cos (^ In (1 + z)) + ^ sin (^ In (1 + z)))xo+ 
(cos (^ In (1 + z)) + ^ sin (^ In (1 + z)))yo) ^ 
b) C = g, yo = aa:o : 

wx{z) = -1 - 1(1 + z)3/2 (^(1 + a) COS (^ In (1 + z)) + ^(5 + a) sin (^ In (1 + z)))xo- 
-|(1 + z)3((l + a) cos (^ In (1 + z)) + ^(5 + a) sin (^ In (1 + z)))^xg 



'weak'; 2.5 < \nB12 < 5 - 'moderate'; lni?i2 > 5 - 'strong'. 

The values of Bayesian evidence for models with wx{z) defined in Table [1] were obtained 



using a nested sampling algorithm [4l||, which iniplementation to the cosmological case 
is available as a part of the CosmoMC code j42l . Si, called CosmoNest 144144 711. It was 
changed for our purpose. We assume flat prior probabilities for the model parameters in 
the following intervals: fim,o ^ [0,1] and Wq G [—2,0], Wi G [—3,3] (Model 1); Wq G [—2,0], 
G [0,2] (Model 2a and Model 3a); Wq G [-2,2], G [0,2] (Model 2b and Model 3b); 
xo G [-1,1], ?/o G [-1,1] (Model 4a); Xq G [-1,1], a G [-3,0] (Model 4b). The values of 
evidence were averaged from the eight runs. 



B. Analysis with SNIa, CMB R and BAO A data 



To compare mode ls g athered in Table [T] we use information coming from the sample of 
iVi = 192 SNIa data (isj, which consists of the ESSENCE sample ^ and a SNIa detected 



by HST [50(]. After suitable calibration SNIa could be treated as standard candles and 
tests on the assumed cosmology could be done. In this case the likelihood function has the 
following form 



£gjyf oc exp 




theor , obs^2 



0"f 



(7) 
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where cxj is known, = rrii — M (mj-apparent magnitude, M-absolute magnitude of 
SNIa), nf"""' = SlogioI^Li + M, M = -Slogioi^o + 25 and = Hodu, where is 
the luminosity distance, which (with the assumption that the Universe is spatially flat) is 
given by du = (1 + Zi)c f^' -^f^ and H{z) is defined in equation (j2]). After an analytical 
marginalisation over the nuisance parameter Ai one can obtain the likelihood function £sn 
which does not depend on the parameter Hq. 

We also include information coming from the CMB data using measurement of the shift 
parameter = 1.70 ± 0.03 for z^ec = 1089) 0, l5l[ |. which is related to the first acoustic 



peak in the temperature power spectrum and is given by = Q J^"^"" -j^dz. Here 

the likelihood function has the following form 



Cfi oc exp 



^^thcor ^obs^2" 



(8) 



As the third observational data we use the SDSS luminous red galaxies measure- 
ment of A parameter (A"^*^ = 0.469 ± 0.017 for za = 0.35) [52|, which is related to 
the baryon acoustic oscillations (BAO) peak and defined in the following way A^^'^"^ = 

_ 1 2 

^/flra,o (^ ^Ho^ ^ ^ 7X lo^ '^)^^ ^ ' ^^^^ '^^^ likehhood function has the following 
form 

^^theor ^obs'j2 



Ca oc exp 



2a\ 



(9) 



The final likelihood function used in analysis is given hy C = Csn'^R'^a- 
The results, i.e. values of In Bu together with their uncertainties, computed with respect to 
the model with linear in a parameterisation of wx{z), are presented in the first column of 
Table [ml 

As we can conclude there is weak evidence to favour model with purely oscillations (2b) 
over the model with linear in a parameterisation. The comparison with models 2a, 3b and 
4a is inconclusive, which means that the data set used in analysis is not enough powerful 
to distinguish those models. Additional information coming from different data set or more 
accurate data set is required. The value of logarithm of the Bayes factor calculated with 
respect to model 4b is close to —1, which could indicates on the weak evidence in favour of 
it, but more information is needed to make the conclusion more robust. There is moderate 
evidence to favour model 1 over the model with damping oscillations (3a). 

We can also check if the damping term, i.e. (1 + 2)^, is required by the data. Comparing 
model 2a with 3a one can conclude that there is moderate evidence to favour the purely 
oscillations (lni?2a3a = 4.32), while comparing 2b with 3b that there is weak evidence to 
favour purely oscillations (In52b3b = 2.13). 

The comparison of the best parameterisation from the set of models with purely oscilla- 
tions (2b) over the best one among the models with damping oscillations (4b) does not give 
conclusive answer: lni?2b4b = 1-06. 

Finally one can compare models with dynamical dark energy with parameterisations 
of the equation of state gathered in Table [T] with the simplest alternative, i.e. the ACDM 
model with w = —1. As one can conclude this model is still the best one in the light of 
SNIa data, CMB R shift parameter and BAO A parameter. However the conclusion from 
the comparison of the ACDM model with the purely oscillating model (2b) is inconclusive 
(ln5ACDM,2b = 1-05). 
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FIG. 3: The functions w{z) for the ACDM model, model with linear in a parameterisation of 'w{z) 
(1), model with purely oscillations (2b) and model with damping oscillations (4b), calculated for 
the best fit values of model parameters (SNIa+CMBR+BAOA data). 

It is interesting that the model with purely oscillations (2b) is favoured by the data 
over the model with linear in a parameterisation of w{z). Also the model with damping 
oscillations (4b) fares well when compared with model 1. One can try to understand the 
reason of such conclusions. In Figure [3] the functions w{z) for the ACDM model, model 
1, model 2b and model 4b are presented, which were calculated for the best fit values of 
the model parameters (obtained in the analysis with the SNIa, CMB R and BAO A data). 
While in Figure H] we present the corresponding distance modulus vs redshift relations (with 
the additional assumption that Hq = 72 kms~^ Mpc~^). 

As we can conclude in spite of the prominent differences in the functions w{z), the 
distance modulus relations are nearly identical for those models. One should keep in mind 
that parameters were fitted using the data which are based on the luminosity distance. In 
this case w{z) is integrated twice. 

C. Analysis with H{z) data added 

It is interesting to consider the relation H{z), as it depends on the w{z) through one 
integral. Unfortunately the present Hubble function measurements on different redshifts are 
small and inaccurate. However H{z) data set could gives us another insight into the problem 
considered. 

One possibility to measure the Hubble parameter as a function of redshift is based on the 
differential ages dt/dz of passively evolving luminous red galaxies (LRG), which correspond 
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FIG. 4: The distance modulus vs redshift relations for the ACDM model, model with linear 
in a parameterisation of w{z) (1), model with purely oscillations (2b) and model with damping 
oscillations (4b), calculated for the best fit values of model parameters (SNIa+CMBR+BAOA 
data) and with the assumption that Hq = 72 kms~^ Mpc~^. The SNIa data set is also presented. 



to the Hubble function through the relation 



H(z) 



1 dz 



1 + z dt 



(10) 



Using Gemini Deep Deep Survey and archival data the authors of [53| obtained nine values 
of the Hubble parameter for different redshifts in the range 0.09 < z < 1.75. Although this 
data set is small and has large uncertainties we include it in our analysis. 

Another method to determine the Hubble function values at various redshifts is based on 
the line of sight (LOS) baryon acoustic oscillation scale measurements. The scale of the BAO 
in the radial direction depends on the H{z). On the other hand the precise measurement of 
this scale is given by the CMB observations, so the comparison gives us the value of Hubble 
parameter. Based on this method and using the SDSS DR6 luminous red galaxies data the 
authors of [s^] obtained the values of H at three different redshifts. The uncertainties are 
highly reduced when compared with the previous data set. We include those points in our 
analysis. 



55 



To complete our H(z) data set we use the HST measurement of Hq 
We repeat previous calculations with the additional iV2 = 13 Hubble function measure- 
ments. The corresponding likelihood function has the following form: £ = Csn'C.rCa'C.h, 
where 



Cff (X exp 



thcor / 



h: 



obs\2 



a 



Hi 



(11) 
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FIG. 5: The H{z) functions for the ACDM model, model with linear in a parameterisation of w^z) 
(1), model with pm'ely oscillations (2b) and model with damping oscillations (4b), calculated for 
the best fit values of model parameters (SNIa+CMBR+BAOA+H data). The filled square, circle 
and filled circle points correspond to observational H data from |53|], 5j] and 55|], respectively. 



The values of In Bu together with uncertainties are gathered in the second column of Ta- 
ble mil As one can conclude the inclusion of H[z) data does not change our conclusion in 
most cases. There is still weak evidence to favour of model with purely oscillations (2b) over 
the model with linear in a parameterisation of w. Evidence in favour of model 4b becomes 
slightly greater (weak evidence to favour this model over the model 1). The evidence against 
model 3a is even greater than in previous calculations, we find strong evidence against it. 

The evidence to favour model with purely oscillations (2a) over the model with damping 
term (3a) is strong {B2a3a = 6.2), while the evidence in favour of model 2b over the model 
3b is moderate (-B2b3fe = 2.57). 

The comparison of the best model among the models with purely oscillations, i.e. 2b, 
with the best model from the set of models with damping term, i.e. 4b, does not give the 
conclusive answer B2bAb = 0.99. As one can conclude the ACDM model is still the best one 
from the models considered, however the conclusion from the comparison with model 2b is 
inconclusive (-BACDM,2fe = 0.97). 

We present the relations H{z) for model 1, model 2b, model 4b and the ACDM model 
in Figure El The Hubble functions were derived for the best fit model parameters in the 
analysis with SNIa, CMB R, BAO A and observational H(z) data. 

As one conclude the relations H{z) for models considered are similar in the redshift range 
under consideration. More data with better quality is required. The most promising future 
H data will come from the BAO measurements. This method gives us much more precise 
data points when compared with the alternative method. 
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TABLE II: The values of distortion parameter /3, bias parameter b and corresponding growth rate 
function f = f3b which are used in the calculations. 



z 


/? 


b 


/ 


references 


0.15 


0.49 ±0.09 


1.04 ±0.11 


0.51 ±0.11 


[52, [58] 


0.35 


0.31 ±0.04 


2.25 ±0.08 


0.70 ±0.18 


[5] 


0.77 


0.70 ±0.26 


1.30 ±0.10 


0.91 ±0.36 


[59] 


3.00 






1.46 ±0.29 


[60] 



D. Analysis with growth rate function data added 



The conclusions stated before are based on the geometrical dark energy probes. It is 
interesting to check how the inclusion of the dynamical probes, related to the growth of 
structures, will change the results. We consider observations of the growth rate function /, 
which is related to the growth function D by the following formula f = dlnD/dlna. Its 
evolution in the general relativity framework is described by the following equation 

= -/'-/ - ^(1 - ^M)w{a)^ + ^nM, (12) 

where Qm(a) 



The values of growth rate at various scale factor (a) for considered models were 

obtained with the help of eq. [T2j It was solved using numerical methods, with the assumption 
that /(a ~ 0) = 1. 

The observational growth rate data {f°^^) could be obtained through the measurements 
of the redshift distortion parameter /3. It is observed through the anisotropic pattern of 
galactic redshifts on cluster scales. It is related to the growth rate function by the following 
formula P = f /b. The so called bias parameter b reflects the fact that the galaxy distribution 
does not perfectly trace the matter distribution in the Universe. Currently there are only 
few measurements of / available (see Table HTj) . This data set is similar to the one presented 



m 



56| . We do not consider the data points aX z = 0.55 and z = 1.4, as the bias parameter 
was derived with the help of the value /3 in those cases. The measurement at z = 3 was 
obtained in different method, which does not rely on /3 and b parameters. The value of / is 
finding in the analysis with Ly-a forest data. 

It should be kept in mind that this data set was obtained with the assumption of the 
ACDM model. Its inclusion in the analysis with other models could decrease its reliability 
and the results should be treated with care. 

The likelihood function used in analysis is of the following form £ = CsN^^R^^A^^Hi^f-, 
where 



Cf (X. exp 



(13) 



where Nj, = 4. The values of In Bn and its uncertainties are gathered in the third column 
of Table [ml 

As one can conclude the final conclusions do not change in all cases. The data set is not 
informative enough to change results. 
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FIG. 6: The f{z) functions for the ACDM model, model with linear in a parameterisation of w{z) 
(1), model with purely oscillations (2b) and model with damping oscillations (4b), calculated for 
the best fit values of model parameters (SNIa+CMBR+BAOA+H+f data). 



TABLE III: The values of ln(Si,;) = \nP{D\Mi)-\nP{D\Mi) calculated with respect to the model 
with linear in a parameterisation of wx{z) for different data sets. 



MODEL 


SNIa+CMBR+BAOA 


SNIa±CMBR±BAOA±H 


SNIa±CMBR±BAOA±H±f 


1 











2a 


0.29 ±0.22 


0.33 ±0.21 


0.18 ±0.22 


2b 


-2.08 ±0.13 


-2.14 ±0.20 


-2.24 ±0.23 


3a 


4.61 ±0.18 


6.53 ±0.24 


6.3 ±0.23 


3b 


0.05 ±0.12 


0.43 ±0.23 


0.34 ±0.25 


4a 


0.5 ±0.13 


0.31 ±0.23 


0.25 ±0.24 


4b 


-1.02 ±0.18 


-1.15 ±0.23 


-1.21 ±0.23 


ACDM 


-3.13 ±0.16 


-3.11 ±0.22 


-3.3 ±0.23 



In Figure [6] one can find a plot of the growth rate as a function of the scale factor for the 
ACDM model, model 1, model 2b and model 4b, calculated for the best fit values of model 
parameters (in the analysis with SNIa,CMB R, BAO A, H and f data). 

The relation /(a) for model 4b differs from the other relations. Anyway the data points 
have large uncertainties, which prevent this set to distinguish models considered. This is in 
agreement with our previous conclusion. 
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E. Discussion 

In spite of the fact that models with oscillating relation for w{z) (i.e. 2b and 4b) fare 
well when compared with the model in which w{a) is a linear function of scale factor, the 
oscillating behaviour is not seen in the w{z) vs z plots in the redshift interval considered. 
While the frequency parameter of the model 4b is fixed by the theory, it appears as a free 
parameter in model 2b. The relation for w{z) in model 4b is complicated. However it can 
be rewritten as sum of sine and cosine components, with the amplitudes, which depend 
on z, as well as on the model parameters. On the other hand frequency parameters Wg 
are equal to \/7/2 or \/7. If we consider the relation |i(7sln(l + z)\ — 2%, we can claim 
that the oscillating behaviour should be observed in the redshift range of about Az ^ 10. 
Unfortunately most of the data points used in analysis are for z < 2. The oscillating 
behaviour is not seen. More observations at higher redshifts are needed. On the contrary, 
as was stated before, the frequency parameter of model 2b is a free one. The assumed prior 
range for this parameter (i.e. Ws G [0,2]) corresponds to the period of oscillation of at least 
Az ~ 22. It is again too big to be observed with the present data sets. It is interesting 
to consider the situation in which the oscillating behaviour could be detected. It can be 
done by assessing a different prior range for frequency parameter. We repeat calculations 
for model 2b, with the assumption that Wg G [2,4.5] (model 2bl). It corresponds to period 
of oscillations of at least Az ~ 3 (the redshift of the most distant data point, of course 
apart from the one at z = 1089). The value of logarithm of the Bayes factor, calculated 
with respect to model 2b, is equal to \nB2b,2bi = 4.1 ± 0.17. This means that the evidence 
against model 2bl is moderate. We can conclude that available data sets prefer a model 
with period of oscillations larger than could be detected nowadays. 

IV. CONCLUSIONS 

We use the Bayesian method of model selection to compare the FRW models with different 
functional forms of dynamical dark energy (different parameterisations of the EoS). We 
examine two categories of parameterisations: a priori assumed and derived from the model 
dynamics. We show that two parameterisations are favoured over most popular linear with 
respect to the scale factor. 

In particular we obtain following results: 

• parameterisation with purely oscillations, i.e. 2b, is the best one among parameteri- 
sation considered in this paper; 

• there is weak evidence to favour this parameterisation over the linear in a parameter- 
isation of EoS (this conclusion is based on the SNIa, CMB R and BAO A data sets 
and does not change after inclusion observational H and / data); 

• data sets used in analysis prefer model 2b in which oscillating behaviour could not be 
detected nowadays 

• comparison of model 4b with the linear in a parameterisation of EoS does not give 
conclusive answer when it is based on SNIa, CMB R and BAO A data, but after the 
inclusion of H and / data we find the weak evidence in favour of this model 
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• the comparison of the ACDM model with the model with dark energy parameterised 
as 2b is inconclusive, a more accurate data set is required to distinguish those models; 

• damping term, i.e. (1 + 2;)^, which appears in parameterisations (3) is not supported 
by the data used in analysis. 

In study of cosmological constraints on the form of dark energy the most popular methodol- 
ogy is study of the viability of different parameterisations for the equation of state parameter. 
They are postulated rather in the a priori forms without connection with true model dy- 
namics. Our approach is different because we claim that if model dynamics is closed then 
corresponding form of dark energy parameterisation should be forced. It is because we tested 
the FRW model with dark energy rather than the parameterisation w{z) itself. 

Appendix A: Linearised formulas for w(z). 

Here we present linearised formulas for w(z) around the critical point corresponding to 
the deSitter state for monotonic and oscillating evolution toward this point [22,]: 

„,on^ -il-SO + m,a)a--'/' + m,a)a-^ 

(1 - 30 + 6^1 - GOM^a'"' + Ba-''^)a-y^ + 3^(1 - 60(^a"' + Ba-'^^fa-^' 

(Al) 

where = = ^^1^, A = \xo + V^^J^§i(\xq + |yo), B = \xo - 

a/3-^ {\^o + Iz/oj, and are the initial conditions for ip and ip', respectively, 
and /i = 2e^o((3(l - 4^ - 4az(l - 30) Aa"' + (3(1 - 40 + 4a/(l - 30)5a-"'), h = 
(-1(3-40 + 1501-20) (^a"'+5a""')Vai(3(l-40-80l-30) (^^a^"' -E^a'^"') - 



af{l - 40 (^Aa"' - Ba'"'^ , and 



osc ^ -(1 - 30 + ^1(0 a)a-^^' + ^2(e, a)a-' . , 21 

X (1 _ 3^) + 6^(1 _ 60^o/i(e, a)a-^/^ + 3^1 - QOh'i^, a)a-^ ' ^ ^ 

where h = xq cos (oosc In a) + ^ sin (oosc In a) (|xo + |yo) , fl'i = '2^ipo(^{l - 6^)h - 
4(1 - 30 ((2:0 + 2/o)cos(aosclna) - aoscS^o sin (oosc In a) - sin (oosc In a)(|a;o + |?/o))) , 
92 = 01 - 60^^ - (1 - 20(1 - 60 (z/o cos (oosc In a) - Oosca^o sin (ttosc In a) - 
2^ sin (aosc In a) (^Xq + |?/o)) ^ - 401 - 30 ({xq + yo) cos (oosc In a) - Oosca^o sin (oosc In a) - 

sin (oosc In a){^XQ + ^2/0) j ? where ttosc = -^y^^i^ and Xq, yo and tpo have their usual 
meaning. 

Note that in all cases purely oscillating scenario does not exist. 
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